Strong-coupling properties of unbalanced Eliashberg superconductors 
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In this paper we investigate the thermodynamical properties of "unbalanced" superconductors, 
namely, systems where the electron-boson coupling A is different in the self-energy and in the Cooper 
channels. This situation is encountered in a variety of situation, as for instance in d-wave supercon- 
ductors. Quite interesting is the case where the pairing in the self-energy is smaller than the one in 
the gap equation. In this case we predict a finite critical value Ac where the superconducting critical 
temperature Tc diverges but the zero temperature gap is still finite. The specific heat, magnetic 
critical field and the penetration depth are also evaluated. 
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I. INTRODUCTION 

The Eliashberg's theory of superconductivity repre- 
sents an elegant and powerful formalism to extend the 
BCS theory to real materials. Main achievements of the 
Eliashberg's theory are the generaUzation of the BCS 
scheme to the strong coupling regime, where the dimen- 
sionless electron-phonon coupling constant A can be of 
the order or larger than unit, A ^ 1, and the inclusion of 
the retarded nature of the electron-phonon interaction, 
characterized by the phonon energy scale ujpi^. 1 Within 
this framework it was possible to understand and predict 
a number of characteristic features of the strong coupling 
regime, as a 2A/Tc ratio larger than the BCS limit, the 
temperature dependence of the magnetic critical field and 
of the specific heat, the appearance of phonon features 
for w > A in the tunneling and in the optical spectra^^ 
Some results of the Eliashberg theory have become widely 
used paradigmatic milestones, as for instance the employ- 
ment of McMillan- like formulas, Tc « exp[— (1-|-A)/A], to 
estimate the critical temperature and its dependence on 
the microscopic interaction in generic superconductors 
The strong-coupling limit Tc, A ^ cjq of the Eliashberg's 
theory has also been examined in details, showing a dras- 
tic change of the superconducting properties, with, for 
instance, Tc, A cx y/XuJo^>^>^ 

Different evolutions of the Eliashberg's theory have 
also been later introduced in the course of years to adapt 
it to the particular cases of specific materials. Multi- 
band effects,— iSiii^i^ii^ anisotropy and non s-wave sym- 
metries of the order parameterj^ d^d'^d^d^d'^ effect of 
vertex correction s ^^d^i^" have been for instance consid- 
ered. In all these cases one should consider in prin- 
ciple the possibility that the electron-phonon coupling 
(or any kind of other mediator) can be substantially 
different in the self-energy and in the superconducting 
Cooper channels. This is most evident in the case of 
d-wave pairing. For example, if we assume a factorized 
interaction, a^F(k, k'; w) = a'^F{u!)ipdO<)'4'dO^'), where 
a^T'(k, k'; uj) is the generic anisotropic Eliashberg's func- 
tion, and V'd(k) is the wave-function for the d-wave sym- 
metry, we would get no contribution in the self-energy. 



This is expected for instance in the case of a spin- 
mediated coupling where the exchange energy J(k, k') 
is factorizable as « J'0(/(k)'!/'d(k')r^'^ the charac- 
teristic energy scale of the pairing, Wgf, is given by the 
spin-fluctuation spectrum. Of course this is an extreme 
limit, and in real systems there will be finite contribu- 
tions in both the self-energy and the Cooper channels, al- 
though in principle arising from different electron-boson 
modes. In any case, there is no reason to expect that 
the electron-phonon coupling Xz relevant for the Z{lj) 
renormalization wave-function to be the same as the one 
Aa ruling the gap equation. 

In this paper we investigate in details the consequence 
of a different coupling in the Z{uj) wave-function and 
in the gap equations. We define this situation as "un- 
balanced" Eliashberg theory. We focus here on ther- 
modynamical quantities which can be evaluated in the 
Matsubara space. Spectral properties, involving analyt- 
ical continuation on the real axis, will be investigated 
in a future publication. We show that, contrary to the 
common feeling, an unbalanced coupling in the Eliash- 
berg's theory has important and drastic differences with 
respect to the conventional Eliashberg phenomenology. 
In particular we show that for Xz < Aa the supercon- 
ducting critical temperature Tc is strongly enhanced for 
finite values of A ~ 1, and in the infinite bandwidth limit 
Tc even diverges. We also show that these new features 
are strictly related to the retarded nature of any boson 
interaction, accounting for the fact that this phenomenol- 
ogy was never discussed in the case of the non-retarded 
BCS theory. 



II. CRITICAL TEMPERATURE VS. A 

Let us start by consider the Eliashberg's equations for 
the simple representative case of an Einstein boson spec- 
trum. To simplify the notations, we define rj — A^/Aa 
the ratio between the electron-boson coupling in the Z 
and in the gap equations, and we simply denote A = Aa- 
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In the Matsubara space we have 



UJ„ 



(2) 



where is the energy scale of a generic bosonic me- 
diator. Eqs. HI)-® can be easily generalized in the 
case of a d-wave symmetry for the gap order parameter 
A„ — > A„'0d(k) in the Cooper pairing. In the weak- 
intermediate regime, defined as Tc/wq, A/cjo ^ 1, a 
simple analytical solution for Tc and A is provided by 
the square-well model. ^ Along this line one recovers, ac- 
cording the common wisdom, a generalized McMillan-like 
formula Tc ~ exp[— (1 + A2)/Aa], which predicts an up- 
per limit for Tc in this case as well as in the perfectly 
balanced r] — \ case. The validity of such result is how- 
ever limited to the weak- intermediate case where Tc/wq; 
A/wq <C 1. In the balanced case, for instance, a careful 
analysis shows that, in the strong coupling regime Tc/loq, 
A/wq S> 1 the critical temperature as well the supercon- 
ducting gap do not saturate for A ^ oo but they scale 
asymptotically as Tc, A oc -/Xujf)^^^^^ 

A first insight that things can be radically different for 
an unbalanced Eliashberg's theory comes from a reexam- 
ination of the strong coupling regime. Plugging Eq. (JT]) 
in ([2]) we obtain for T^. 



AttTc^ 



A„ - 77(u;„/w„)A„ 



(3) 



For rj = 1 the term n — m vanishes in Eq. ([3]), so 
that, for Tc 3> loq, the first contribution in the bo- 



son propagator comes from ijJq/[ujI + {uj„ 



V 



[u)q / An'^T^ {n — m)"^]^ This is no more the case for the 
unbalanced case where the leading contribution comes 
from luq/[uj^ + (ti>„ — cJm)^] « Sn,rn- Eq- & rcads thus: 



1 = 



\2n+l\ ' 



(4) 



Note that the temperature T does not appear anymore 
in Eq. For 77 < 1 Eq. ^ implies that there is an 

upper limit A™'*'' ^ 1/(1 — ry) above which the system 
is always unstable at any temperature with respect to 
the superconducting pairing. A detailed analysis (see 
Appendix [A| shows: 
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AAc 



Ac - A' 



(5) 



where Ac = 1/(1 — ?/)■ On the other hand, for rj > 1, Eq. 
(HI) is never fulfilled signalizing that the limit Tc ^ luq is 
unphysical and Tc must saturate for A — > 00 limit. 
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FIG. 1: Critical temperature as function of the pairing cou- 
pling A for different values of 77. Dashed lines: 77 < 1; dot- 
ted lines: 77 > 1; the solid line is the conventional Eliash- 
berg result for 77 = 1. From upper to lower line: 77 = 
0,0.5,0.8,0.9,0.95,1,1.2,1.5,2. 



We would like to point out that the analytical diver- 
gence Tc ^ 00 for A Ac is strictly related to the infi- 
nite bandwidth model employed in Eqs. ([IJ-lH]). On the 
other hand, in physical systems the presence of a finite 
bandwidth W determines an additional regime, Tc ^ IF, 
where the analytical divergence of Tc at Ac is removed 
and Tc oc AIF^ (Appendix [X)) . In this respect the band- 
width W defines an upper limiting regime for Tc. Since 
in physical systems, however, W is some orders of mag- 
nitude bigger than the bosonic energy scale cjq, in the 
following, for sake of simplicity, we shall concentrate on 
the infinite bandwidth limit W 3> Tc, A,a;o, keeping in 
mind however that the analytical divergences found in 
this case will be removed when finite bandwidth effects 
are included in the regime Tc ^ IF. 

In Fig. [T] we show the critical temperature Tc as func- 
tion of the electron-boson interaction A obtained from 
the numerical solution of Eqs. ([I])-©. We see that the 
conventional Eliashberg case 77 = 1, where Tc scales as 
\/A, represents rather an exception than the rule: for 
77 < 1 Tc diverges at finite values of A determining, for 
each 77, a upper value of A above which the system is su- 
perconducting at any temperature, while, for 77 > 1, Tc 
saturates for A ^ 00 at some value which also is depen- 
dent on 77. We can estimate in this case (see Appendix 
|A|) an upper limit for Tc- 
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LUq 



(6) 



Before to proceed on, we would briefly comment on the 
Tc divergence at finite A for 7; < 1. This result seems to 
contradict apparently the BCS scenario Tc cx exp(— 1/A) 
which can be thought as the extremely unbalanced 77 = 
case where the one-particle renormalization processes are 
disregarded. However, a closer look at Eq. ^ shows that 
a fundamental role in deriving Eqs. Q-® is played 
by the proper treatment of the retarded nature of the 
electron-boson interaction, which gives rise to the cor- 
relation between w„ and ujm within the energy window 
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FIG. 2; Matsubara superconducting gap Am (upper panel) 
and ratio 2 Am /Tc (lower panel) as function of the pairing 
coupling A for the same values of rj as in Fig. [1] 



ojQ. In this sense, neglecting the Z- function for 77 = 
in Eqs. ([I])- ([21) corresponds to a "retarded BCS" the- 
ory. This is quite different from the usual conventional 
BCS framework where the interaction is supposed to be 
non-retarded and the frequencies ojn and are uncor- 
related. This scenario can be achieved in the retarded 
BCS context only in the limit coq —^ 00, which enforces 
the Tc/ujq — !• limit, namely, the weak-coupling regime. 



III. SUPERCONDUCTING GAP Am VS. A 

Interesting enough, the A vs. A behavior in an unbal- 
anced superconductor can be quite different from the Tc 
vs. A. In Fig. [21 we plot the Matsubara superconducting 
gap Am, defined as Am = ImiT^o ^n=Q, and the ratio 
2Am/T'c as function of the electron-boson coupling A. We 
remind that, although Am underestimates the physical 
gap edge obtained by the analytical continuation on the 
real axis, the analytical dependence of these two quanti- 
ties is usually the same, so that Am can be reasonable 
employed to study the limiting behavior of the supercon- 
ducting gap in the strong-coupling regime. Detailed in- 
vestigations on the real axis are however needed to assess 
in a more formal way this issue. Fig. [5]shows that, while 
for r] > 1 Am has a saturating behavior similar as Tc, 
in the case 77 < 1 the superconducting gap does not di- 



verge at some finite value of A, as Tc, but rather increases 
linearly with the electron-boson coupling. This different 
behavior can be also understood by applying some an- 
alytical derivations properly generalized for unbalanced 
superconductors.^^ Assumiiig Am 3> luqj and following 
Eqs. (4.29)-(4.35) of Ref. [j, the superconducting gap 
Am is determined by the following relation: 



1 + r;A 



-n-LUo 



ClTjX- 
M ^M 



2A 



M 



-2A|^, (7) 



where ci, C2 are constant factors whose value is dis- 
cussed in Appendix [XI For rj = 1 the terms oc wq/Am 
on both the left and right sides cancel out, so that 
Am oc •v/Awq'^'^ This is no longer true for 7/ 7^ 1. In 
particular, for 77 < 1 we find 



A(l-ry) 



TTUJo 



(8) 



while, for 77 > 1, Eq. ^ does not admit solution signal- 
izing, once more, that the initial assumption Am 3> luo is 
inconsistent in this limit and that Am must be saturate 
for A — !■ cx). By taking into account higher order terms 
we found an upper limit Am for 77 > 1 in the regime 
A ^ 00 in similar way as done for Tf.: 



iM.i 



4(7?Ci - C2)UJQ 

27r(77-l) 



(9) 



Note that Am, max in Eq. Q diverges as 1/ (77— 1) whereas 



Tcmax in Eq. ^ scales as Te n 



l/V^FT. This 



means that the ratio 2Am/Tc is not bounded for 7/ > 1 
and it can be even larger than in the Eliashberg case 
7/ = 1 and formally diverging for 7; ^ 1, in agreement 
with the numerical results shown in Fig. [2l 



IV. TEMPERATURE DEPENDENCE OF AM(r) 

In the previous sections we have studied the strong 
coupling behaviors of the critical temperature Tc and of 
the zero temperature Matsubara gap Am in the limit 
A ^ 1. We have seen for instance that, in the 77 < 1 case, 
Tc diverges at some critical value Ac, so that for A > Ac 
the system is superconducting at any temperature. This 
behavior is in contrast with the one of the Matsubara gap 
Am which is always finite for any A and scales linearly 
with A for 77 < 1. As shown in Fig. [21 these different be- 
haviors are reflected in a ratio 2Am/Tc smaller than the 
BCS limit 3.53 and vanishing for A —^ Ac- In this situa- 
tion an interesting issue to investigate is the temperature 
dependence of the superconducting gap A(T), which is 
reflected in a number of observable physical behaviors, 
as the temperature profile of the magnetic field ffc(T), 
of the London penetration depth Al(T) or of the specific 
heat Cy(T). Also intriguing is the situation 77 < 1 and 
A > Ac, where a finite superconducting gap exists at zero 
temperature but where no finite critical temperature is 
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FIG. 3; Temperature evolution of the Matsubara supercon- 
ducting gap AM(r)/AM(0) for different unbalanced cases. 
Left panel: rj — (Ac = 1) and different coupling (from top to 
bottom) A = 0.5,0.8,0.9,0.95,0.98,0.99,0.995. Right panel; 
rj = 2 and (from bottom to top, but barely distinguishable) 
A = 5,10,20. Inset: same quantities as in the left panel 
{rj — 0) but as function of T/luq. 



predicted. In this case the temperature behavior of the 
gap itself is not clear and needs to be investigated. 

In Fig. [3] we show the temperature dependence of 
Am(T)/Am(0) (defined as Am(T) = A„=o) for differ- 
ent characteristic cases, namely for rj = and different 
A < Ac, and for rj = 2 and different A > 1. Most reg- 
ular is the Am{T) vs. T dependence for > 1, where 
Am(T') follows a conventional behavior, independently of 
the coupling A. This regular behavior can be understood 
by reminding that for 77 > 1, even for very large cou- 
pling A, the values of and of the superconducting gap 
Am (7") are always finite and (at most) of the same order 
of the energy wq. Quite different is the case of 77 < 1, 
here represented by 77 = 0, where Am(T') shows a tem- 
perature dependence remarkably different from the BCS 
one. For A close to Ac, in particular, the superconducting 
Matsubara gap has a first initial drop followed by a more 
regular dependence. This change of curvature represents 
the crossover between a small temperature {T/ujq ^ 1/4) 
to a large temperature {T/uq 1/4) regime, as shown in 
the inset of Fig. [3] where we plot Am(T')/Am(0) as func- 
tion of T/ojQ. Note that, while the value of the critical 
temperature Tc is strongly dependent on the coupling A, 
the initial dependence of Am(T) is only weakly depen- 
dent on A. We remind indeed that the Tc divergence for 
A — > Ac is essentially a by-product of having reached the 
T/ujQ ^ 1 for a finite A in the unbalanced case. We 
can thus understand the results of Fig. [3] in the follow- 
ing way: for low temperature {T/loq ^ 1/4) the super- 
conducting gap probes a pairing kernel which is actually 
increased by the lack of unbalance, but still regular, (re- 
mind that Am(T — 0), contrary to Tc, does not diverge 
at Ac, but it steadily increases as oc A in the strong- 
couphng regime). For low temperature Am (7")/ Am (0) 
follows thus a standard-like behavior which would close 
to some finite not diverging at Ac. When A is however 



close enough to Ac, the range T/loq ^ 1/4 is achieved 
before is actually reached; in this regime high tem- 
perature effects become dominant in the pairing kernel, 
reflected in a change of the AmC^) vs. T trend and in a 
final, physical, Tc which diverges at X —^ Xc- 

It is interesting to investigate also how the super- 
conducting gap Am{T) closes at Tc. In the conven- 
tional, perfectly balanced, Eliashberg theory the normal- 
ized gap A{T)/ttTc (or equivalently A(T)/A(T = 0)) 
scales indeed for T ^ Tc as {A{T) / ttTc)'^ ~ cS, where 
6 — I — T/Tc and where c is a finite constant which, in 
the weak-coupling BCS limit A ^ 1, is c = 0.95 whereas 
for A ^ 1 one gets c = 2. This scenario is qualitatively 
different in the case of unbalanced 77 < 1 superconduc- 
tors (Fig. O left panel) where the constant c strongly 
depends on the coupling A. 

For 7/ < 1 a first insight about the temperature depen- 
dence of Am (T) close at Tc is gained simply by consider- 
ing that, for A — > Ac, Am(0) is finite while Tc 00, with 
a ratio 2Am(0)/Tc 0. This means that, as A ^ Ac, 
the constant c must vanish. This result can be shown in 
a more quantitative way (see Appendix [X| by employing 
a one-Matsubara-gap approximation^ which, for T ~ Tc 
and for Tc 3> wqj is quite reasonably justified. For generic 
77 one obtains thus 



7r2T2 



4(5 



1- A(l-77) 

2- A(l-77)' 



(10) 



where S = l-T/Tc. Note that Eq. HU]) reduces to the 
standard relation (Aq/ttTc)^ = 25 for 77 = l.'^ On the 
other hand, for 77 < 1 



7r2T2 



= 4<5 



Ac -A 
2Ac- A' 



showing that the coefficient c vanishes as c cx Ac 
A^ Ac. 



(11) 
A for 



V. OTHER THERMODYNAMICAL 
QUANTITIES 

The anomalous temperature dependence of Am(T) 
is reflected also in other thermodynamical, measurable 
quantities, as the specific heat Cy(T) or the magnetic 
critical field Hc{T). In order to investigate these prop- 
erties we evaluate numerically the free energy difference 
AT between the superconducting and the normal state^^ 



AT(T) = -7rTA^(0)' 



(12) 



where iV(0) is the electron density of states at the Fermi 
level and Zf^, are the Z-renormalization functions 
calculated respectively in the superconducting and in the 
normal state. From Eq. p2p we obtain the magnetic 
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derivation of ACv, we obtain thus 



FIG. 4: Specific heat Cv{T) for tlie two cfiaracteristic unbal- 
anced cases previous considered. Left panel: rj = (Ac = 1) 
and A = 0.5, 0.8, 0.9, 0.95, 0.98, 0.99, 0.995. Right panel; r] = 2 
and A = 5, 10, 20. Inset: specific heat for = as function of 
T/luq. 



critical field Hc{T) = ^/^8ttAF(T) and the specific heat 
difference as ACy(r) = -Td^AF/dT^. One can also 
obtain the total specific heat Cy (T) ACy (T) + by 
adding the contribution that the system would have in 
the normal state, = -fT, where 7 = (2/3)7r2A^(0)(l + 
77A) is the Sommerfeld constant. In similar way we can 
evaluate the London penetration depth Al(T) asi^ 



2ttT- 



A? 



e^vlN{0)\l{T) 



Zn [col + Al] 



213/2 



, (13) 



where e is the electron charge and vp the Fermi velocity. 

In Fig. 0] we show the temperature dependence of the 
specific heat as function of A for the two representative 
cases ?7 = 0, 77 = 2. In this latter case the specific heat has 
a regular activated behavior and it is almost independent 
of A, in agreement with the corresponding weakly A de- 
pendence of the temperature behavior of Am (T) reported 
in Fig. [31 Note also that for this value of 77 the asymptotic 
value A — > 00 of the specific heat jump ACy /jTc ~ 4.2, is 
quite larger than the BCS hmit, (ACv/7Tc)bcs - 1-43, 
pointing out that the superconductors is in an effective 
strong-coupling limit. We should remark that the asymp- 
totic value limx^oo ACv /"jTc is actually dependent on 
the specific value of the parameter 77 > 1. 

Quite anomalous is also the temperature behavior for 
rj < 1. In this case we see that approaching A — > Ac 
the jump is remarkably reduced. Eqs. (fTU)) - PT|l can be 
used to estimate the jump ACy at Tc of the specific heat. 
Using once more the one-Matsubara gap approximation, 
and employing a standard analysis;^ one can show that 
the formal expression for the free energy difference AF 
close to Tc is just equal as in the standard case. 



AF 



iV(0)(7rrc 



47r2T2 



Ao(r) 



(14) 



ACy 



24 



1 + ??A 



1 + r/- 



A. 



(15) 



Plugging Eq. (fTU)) and ^ in (fTil) , and using the standard 



showing that, for 77 < 1, also ACvj^T^ as (Ao/TrTc)^ 
scales as oc Ac — A for A — > Ac. Once more, Eq. ([TS]) re- 
duces to the standard analytical result ACy /^Tc — 12/A 
for 77 = 1 and A 00. It is interesting to note 
that the vanishing of ACv llT^ for A ^ Ac is mainly 
due to the vanishing of ACy oc (Ac — A)^''^^ whereas 
Tc cx 1/VAc — A. This means that, contrary to the 
perfectly balanced case rj — \ where the vanishing of 
ACyj^Tc is driven by 7rc ~* 00 which overcomes the 
divergence of ACy, in the unbalanced 77 < 1 case the 
specific heat jump ACy is itself vanishing. Such obser- 
vation points out that, although Tc is much higher, the 
superconducting properties in the T/wq regime of 77 < 1 
unbalanced superconductors are much weaker than the 
usual. 

This scenario is once more outlined in Fig. 21 where 
we see that the vanishing of ACy is accompanied by the 
developing of a shoulder at T sa wo/4 (see inset). Above 
this temperature the specific heat scales roughly linearly 
with T as a normal ungapped metal with a smaller and 
smaller jump at Tc as A Ac. 

A similar trend is observed in the study of the temper- 
ature dependence of the magnetic critical field and of 
the London penetration depth Al, shown in Fig. |5l For 
77 = 2 we find again that the temperature dependence 
of both H,\T) and Al^(T) presents again a conventional 
behavior with even a more marked curvature with respect 
to the BCS curve. This is compatible with the specific 
heat jump which is also larger than the BCS limit. Quite 
different is the case for 77 = 0. Here we observe a change 
of curvature at low temperature which is more marked as 
A — > Ac.^^ Once again such change of curvature occurs for 
T'^Lo^jA: and for A ~ Ac it is refiected in a sudden drop 
of ffc(T) and Al^(T) signalling once more that, although 
the critical temperature is strongly increasing, the high 
temperature superconducting properties of unbalanced 
systems are quite poor. 



VI. SUMMARY AND DISCUSSION 

In this paper we have investigated the properties of 
"unbalanced" retarded superconductors, namely systems 
where the electron-boson coupling \z in the one-particle 
renormalization function is different than the one rele- 
vant for the gap equation, Aa- We have shown that 
the superconducting properties are strongly dependent 
on the ratio 77 — Xz/Xa- We have analyzed both the 
cases 77 > 1 and 77 < 1. In the first case we show that, 
contrary to the perfectly balanced case, the critical tem- 
perature is always finite and it even saturates at a finite 
value for A — s- 00. The superconducting properties in 
this case are quite similar to the conventional case in 
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FIG. 5: Temperature dependence of the magnetic critical 
field Hc{T)/Hc{0) and of the London penetration depth 
/ ^l{T) for the two characteristic unbalanced cases pre- 
vious considered. Left panels: rj = (Ac — 1) and A = 
0.5, 0.8, 0.9, 0.95, 0.98, 0.99, 0.995. Right panels; 77 = 2 and 
A = 5, 10, 20. Insets: the same quantities as function of T/uq. 



the weak-intermediate regime where the magnitude of rj 
rules the relevance of the strong-couphng effects. Quite 
anomalous is on the other hand the case 77 < 1, which is 
relevant in a variety of situations as for instance for d- 
wave superconductors. In this case we find that the criti- 
cal temperature, in the infinite bandwidth limit, diverges 
at a finite value Ac « 1 (Ac = 1 for rj = 0). We show 
also that, although for A > Ac the critical temperature 
Tc is strongly enhanced for A — > Ac, the zero temperature 
gap is still finite and it scales as A c>c A. The anomalous 
temperature dependence of the superconducting gap is 
reflected in a variety of other physical properties, as the 
magnetic critical field, the penetration depth and the spe- 
cific heat. All these quantities show a strong anomaly at 
T ^ wo/4 above which the system presents very weak 
superconducting properties, which however persist up to 
the higher critical temperature Tc. We would like to re- 
mark once more that in real systems finite bandwidth 
effects remove the analytical divergence of Tc at Ac when 
T,^W. 

These results suggest an interesting pseudogap sce- 
nario. In the T'^ T < Tc regime indeed, since the 
superconducting binding energies are quite small, the ro- 
bustness of the long-range order of this phase towards 



phase fluctuations, disorder, phase separation and other 
different kinds of instability is highly questionable. In 
case of loss of the long-range order, these weak supercon- 
ducting properties will present themselves as a pseudogap 
phase in the range temperature T'^ T Tc^ where the 
temperature T'^ will appear as the therniodynamical crit- 
ical temperature where true lost range order is lost, and 
Tc ~ Tc will set the pseudogap temperature T*. Note 
that, even in this case, the critical temperature of the 
long-range order will result T^ ^ wo/4, which can be 
significantly higher than the predictions of the standard 
Eliashberg's theory. 

Before concluding, we would like to spend some more 
words about the physical meaning of the Tc divergence for 
•q < 1. On the mathematical ground, we have seen that, 
contrary to the conventional balanced case, in unbal- 
anced superconductors the Cooper instability is driven 
by the n ~ m term in Eq. ([3]). In physical terms this 
corresponds to consider the classical limit T ^ of the 
bosons. This is quite different from the usual Eliashberg's 
theory where only virtual bosons, which a characteris- 
tic energy scale cjq, are responsible of the pairing. Just 
as for the linear behavior p(T) cx AT of the resistivity, 
in the classical limit T ^ ojq the energy scale ujq does 
not provide anymore any upper limit and the effective 
pairing is mainly ruled by the increasing bosonic popula- 
tion n{T) cx T with temperature. In the absence of any 
other competing effect, increasing temperature will result 
thus in a stronger pairing with a positive feedback which 
would lead a high critical temperature Tc ^ wq, and the 
only limiting energy scale in this case is provided by the 
electronic bandwidth. In this scenario a competing ef- 
fect in balanced superconductor is provided by the fact 
that the increase of the boson population, as T increases, 
would act in a similar way in the self-energy channel. 
The competition between these two effects give rise to 
the well-known Tc^/Xujo dependence in perfectly balanced 
superconductors with t] — 1 Such equilibrium does not 
occur however in unbalanced superconductors with rj < I 
where the gain in the Cooper channel is larger than the 
loss in the self-energy. In this case, at sufficiently large T, 
the increase of the Cooper pairing due to the boson pop- 
ulation will prevail over the one-particle renormalization 
effects and a superconducting ordering can be sustained 
up to high temperatures Tc limited only by the electronic 
bandwidth energy scale. 
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APPENDIX A: ANALYTICAL FORMULAS FOR 
ASYMPTOTIC BEHAVIORS Tc»ujo, Am > loo 

In this appendix we derive some useful limit expres- 
sions for the critical temperature Tc and the Matsub- 
ara superconducting gap Am in the regimes Tc ^ loq, 
Am > Wo- 



Critical temperature Tc 



Let us now consider finite bandwidth effects in the rj < 
1 case. For simplicity we shall focus on the representative 
case f] = 0. Finite bandwidth effects can be included in 
the linearized Eliasberg's equations, by replacing the 
term 7r/|wm| —> (2/a;„j) arctan(W^/2cj,„), where W/2 is 
the half-bandwidth in a symmetric particle-hole system. 
In this case, in order to obtain an analytical expression 
for Tc in the limit Tc ^ wq, it is sufficient to retain only 
the n = m term in Eq. ([3]), and, in the one-Matsubara- 
gap approximation, we obtain 



Let us start by considering Eq. ([3]) and let us as- 
sume Tc 3> loq- In this regime we can employ the 
one-Matsubara-gap approximation,^ where the only not 
vanishing terms of the Matsubara gap function A„ are 
A„^o = A„^_i. We have thus: 



AAo -I- AA„i 



-At^Ao^ 



47r2T2 



UJm) 



2Sgn[w„], (Al) 



An = AAn — arctan 



W 
2^c 



(A7) 



Inverting Eq. (|A7)) we obtain thus, for Tc > W^, Tc = 
(M^/27r)/tan(7r/2A) « XW/n^, showing that Eq. (|^ 
is valid as far as ujo Tc W , whereas for Tc 3> 
the divergence at Ac is removed and a linear behavior as 
function of A is achieved. 



2. Superconducting Matsubara gap An 



and, reminding that A_i — Aq and J2m ^guf'^nil^o/I^o" 
(t^„=o - t^m)] = 1, we have 



1 = A(l -r;) + A 



47r2T2 • 



(A2) 



For r] = 1 the first term on the right hand side is zero, 
and we recover the usual result T^=^ = y/\uJo/{2n). On 
the other hand, for 77 < 1 we obtain 



rprj<l 

^ c 

or, equivalently. 



LOq 



A 



27r y 1 - A(l - 77) ' 



rpr]<L _ 



AcA 



27r V Ac - A' 



(A3) 



(A4) 



where Ac = 1/(1 — 77). 

Finally, The same expression (|A3p can be used to eval- 
uated an upper limit for Tc in the case 77 > 1. We can 
write thus 



ujo_ / A 
27rY 1 + A(77- 1)' 



and, for A — > cx). 



rpri>l 
c.max 



Wo 



27rV^^ 



(A5) 



(A6) 



Note that Eq. (|A6[) predicts a saturating upper limit for 
Tc''^^ for A 00, violating the assumption Tc ^ ujq. 
This value, T^JJiax: can be thus assumed as a upper limit 
for the asymptotic behavior of TJ}^"^ in the regime A 
00. 



Let us now investigate the asymptotic behavior of the 
zero temperature Matsubara gap Am. Once more, we 
assume the limit Am ^ ^^o s-nd we shall check later the 
consistency of this ansatz. To get an analytical expres- 
sion for Am, we essentially follow Refs. I3II23I. Transform- 
ing, in the zero temperature limit, the Matsubara sum in 
an integral, Tj^m ~^ /c?w/27r, we can write: 



Z{uj)A{lj) 



Xrj 
2^ 



A 



cLu'lu' 



^Cj'2 + A2(tj')^0 + (^-^') 



dLu'A{uj') 



Wo 



We also employ the simple modelj^i^^ 

w'i^A^' (AlO) 



which we have tested numerically to be appropriate even 
for unbalanced superconductors. Here a is a constant of 
the order of unity. Using this model and expanding Eqs. 
[9]) in powers of wq/Am we end up with: 



1 + r/A 



TTUJQ 



2A 



M 



C177A 



A2 



A 



TTLOo 



2A 



M 



C2A-5-(All) 



where ci, C2 are constant factors depending on a. For 
a = 1 we have for instance ci — 2(\/2— 1/2), C2 — v^,— 
while for a = 2 we have ci = -/S — 1, C2 = -^75/24 

In the perfectly balanced case rj ~ 1 the linear terms 
in Wo/ Am in Eq. (|A11[) cancel out, and one recovers 
the usual result A^ = ^/X^/c^ 

Ci 



C2LO0 (we remind that 
C2 is a positive quantity). On the other hand, for 
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r] < 1 the solution of Eq. (jAlip is dominated by the 
hncar terms giving: 



(A12) 



More complex is the case 77 > 1, where retaining only 
the linear term is not sufficient. In this case one needs 
to consider explicitly also the (wq/Am)^ and to solve the 
quadratic equation (jAltp . We obtain: 



AmV ^ A7r(^_l) Mm 

LOo J 2 \ UJQ 



- A(77ci - C2) = 0(Af 3) 



so that 



A'?>1 
M,max 



A7r(77- f) 



A27r2(?y- f)2 



4A(?7Ci 



C2 



4(7701 - C2) 

27r(r/- f) ' 



(A14) 



showing that for 77 > 1 Am saturates at a finite value for 
A — > 00, and Eq. (IA14P can be considered an upper limit 
for it. 



3. Temperature dependence of Am{T) 

In this section we investigate the temperature behavior 
of the Matsubara gap Am{T) close to Tc. We assume 
once more th limit Tc 3> ujq and thus the validity of the 
one-Matsubara-gap model. Within these approximation 
we can rewrite Eqs. ([T]) as: 



1+77A J2 



2 , / T^sgn w„ 

UJn + (LOn — ^m) 



= 1 



^0 



1 



where we have used Aq = A_i, and in similar way. 

Ao 



ZqAq — XnT- 



^n^T^ + A2 V ^ ■ 47r2T2 



1 



(A16) 



Plugging (jA15|) in (jA16|) we have: 
ttTAo 



Ao = A(l-r;) 



+A(l+77) 



TTTAn 



^712^2 + A2 47r2T2 47r2T 



which we can rewrite as: 



0f2T2 



A(l - 77) + A(l + 77) 



.(A18) 



Expanding left- and right-hand sides of Eq. (|A18p at 
the second order in A and at the linear order oi 5 = 
1 - T/Tc, we have: 



A2 

27r2T2 



2,5 Aw^ 



47r2r2 + j^XojI 



25 



1- A(l-7;) 

2- A(l-77)' 



(A19) 



where we made use of Eq. (jA2[) . Eq. (IA19P reduces to 
the standard relation [Aq/ttTcY — 2(5 in the perfectly 
balanced case 77 = 1.— On the other hand, for 77 < 1 we 
have 



A2 
7r2T2 



45 



Ac -A 
2A,- A' 



(A20) 



which shows that A(r) ~ c^/Tc - T (for T ~ T^) with a 
vanishing coefficient c for A Ac. Finally, for 77 > 1, Eq. 
(|A19[) is well-behaved in the limit A ^ 00 and it gives 



A2 
rr2r2 



4,5, 



(A21) 



with a coefficient c twice larger than the usual. 

Eq. (jAlSp can be employed also to study the temper- 
ature behavior of the superconducting gap in the regime 
77 < 1 and A > Ac, where the system is superconducting 
even at high temperature and no finite is predicted. 
Performing the limit T 3> we obtain: 



A^T2 + A2 



A(l-r;), 



(A22) 



and 



Ao ~ 7rTV(A/Ac)2 ~ 1, 



(A23) 



showing that Ao increases linearly with T in this regime. 
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